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^toract-Haemodialysis has been considered an effective therapy 
for patients with end-state renal disease; however, patients can 
suffer from adverse side effects during haemodialysis. About 
25% of patients who need haemodialysis treatment experience 
acute complications (haemodynamic instability) during the 
treatment sessions, which provoke discomfort in the patients due 
to hypotension, fainting, vomiting, etc., and usually require 
relatively long recovery periods. These instabilities are related to 
the excess shift of fluids between the extracellular (ECW) and 
intracellular (ICW) spaces as shown in several publications 
during the past 10 years. In this work we propose a non-invasive 
method based on local multifrequency bioimpedance 
measurements that allow us to determine the fluid distribution 
and variations during haemodialysis. For that purpose a 
measurement system was developed for the measurement of the 
body fluid balance changes (ratio ECW/ICW) during 
haemodialysis sessions. Clinical measurements were done using 
10 HD patients during 60 HD sessions. Bioimpedance data, 
arterial blood pressure, blood volume and blood heamatocrit 
variations were recorded continuously during the HD sessions. 
Impedance values at infinite and zero (Roo,Ro) frequencies were 
extrapolated from Cole-Cole mathematical model. These values 
are assumed to represent the impedance of total tissue fluid and 
the impedance of the extracellular space respectively. 
Estimators for the extracellular and intracellular fluid volumes 
were developed using Hanai theory of mixtures with values of 
(Roo,Ro)- Significant decrease in the ECW volumes were 
recorded during HD session for all patients, however ICW 
variations were not similar in all patients. In this phase of the 
research work, the selected patients were stable; however, in 
future work non-stable patients will be selected and the 
electrical impedance measurements of the fluid variations will 
be used to develop an indicator of hypovolemic crisis during HD 
thus allowing a better control of the adverse effects during 
dialysis sessions. 

Keywords - Multifrequency Bioimpedance, Haemodialysis, body 
fluid compartments. 

I. Introduction 
Numbers of patients suffering from Chronic Kidney Disease 
(CKD) are very serious. 20 million Americans (1 in 9 US 
adults) have CKD and another 20 million more are at 
increased risk [National Kidney Foundation, 2006]. The 
progress of kidney disease may lead to kidney failure, which 
requires dialysis or a kidney transplant to maintain life. 
Haemodialysis has been considered an effective therapy for 
patients with end- state renal disease. CKD patients usually 
suffer from excess of toxic end-products of nitrogen 
metabolism (urea, creatinine, uric acids, and etc.) 
accumulated in the blood and tissues. Furthermore, these 
patients have a defected salt-water balance that will result in 
shrunken body cell mass with relatively expanded ECW 
space [Jaeger and Ravindra, 1999]. About 25% of patients 
who need haemodialysis treatment can suffer from adverse 



side effects during haemodialysis such as hypotension, 
headache and vomiting. Most of these irritating symptoms are 
related to the excess shift of fluids between the extracellular- 
ECW- and intracellular spaces - ICW- as shown in several 
publications during the past 10 years. The applications of 
bioimpedance techniques in the analysis of body 
compartments are broad. Hundreds of publications in this 
field are available in the literature from the last 30 years. 
Bioimpedance measurements used to evaluate the correlation 
between the Total Body water (TBW) and electrical 
impedance were first reported by Hoffer(1969) using a single 
frequency impedance analysis (BIA) [7] in which he 
presented the correlation between the total body water(TBW) 
and total body impedance. Bioimpedance analysis (BIA) was 
proposed as a method to measure the fat-free mass (FFM) by 
Lukaski (1985) [8]. The applications of BIA were further 
extended to measure the extracellular water (ECW) [9] and 
body cell mass (BCM) [10]. Most Bioimpedance 
measurements are performed using four electrodes: The 
current (<lmA) is passed through the 2- injection electrodes, 
while the voltage drop is measured using the 2-detection 
electrodes. From the injected current and the measured 
voltage drop values, the impedance is derived and the total 
volume can be estimated. Bioimpedance spectroscopy (BIS) 
combines suspension theories for volume estimation [11] 
with biophysical model in which a spectrum of frequencies is 
used to describe resistance and reactance of the whole tissue. 
The resistance of intracellular (ICW) and extracellular ECW) 
mediums can be determined and their volumes are estimated 
using BIS with a range of frequencies (5kHz to IMHz). 

II. THEORY 
A. Bioimpedance as a function of frequency 
In biological systems cell membranes act like capacitors, thus 
the injected current at low frequency (LF) will only penetrate 
the ECW; while at very high frequencies(HF) the current will 
penetrate both the ECW and ICW. 
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Fig. 1. Electrical current paths in biological tissue at LF and HF. 
A useful method for illustrating the behavior of tissue 
impedance as a function of frequency is done by fitting the 
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measured bioimpedance with the Cole-Cole mathematical 
model developed by (Cole and Cole, 1941) [9]. 
The Cole-Cole model obeys the following equation: 

z{w) = R^^ ^;~\_^ (1) 



\ + {jQ)T^) 



Where: 



Z(cd): impedance as a function of the angular 

frequency (radians per second) 

Rq: resistance at co= 

R oo '.resistance at co = oo 

T :time constant 
a (0 <a< 1) is a characteristic parameter of the distribution of 
the relaxation frequencies of the various structures making up 
the heterogonous material. 

As shown in Eq. 1 bioimpedance is complex, therefore it can 
generally be described as : 
Complex: Z = R -\- jX 



Module: Z 



= ^/^ 



+ X' 



(2) 



Phase: 



tan 



"(— ) 
R 



R and X are the real and imaginary part of the complex 
bioimpedance. Bioimpedance module, phase, real and 
imaginary part are usually plotted as a function of frequency 
domain. Fig.2-4 display the behaviour of the impedance in 
several forms, at the frequencies of 10 kHz, 20 kHz, 50 kHz, 
100 kHz, 220 kHz and 500 kHz.(real data measured from 
human muscular tissue). Fig.2 shows the real part of the 
measured data (small circle points) plotted with the fitting 
Cole-Cole model (solid line) as a function of the frequency. It 
is important to note that to measure values close to the 
limiting ones the frequency range would have to be extended 
at least in two decades. 
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Fig. 2. Real part of bioimpedance as function of frequency. 

Fig.3 displays the imaginary part in absolute value because it 
is negative for biological media 
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Fig. S.Imaginery part of bioimpedance as function of frequency. 



The phase angle plotted in Fig.4 displays a behaviour similar 
to the imaginary part. Phase angles are also negative, and in 
this case quite low. To present to Cole-Cole data in the 
complex plane, the real component R is plotted versus 
imaginary component. Fig. 5 displays a complex plane 
diagram, sometimes called the Wessel diagram. A circular arc 
is obtained. The circular arc corresponds to a half circle only 
ifa=0. 
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Fig. 4. Phase angel as function of frequency. 
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Fig. 5. Bioimpedance in the complex plane. 

The ECW resistance (Re) is equal to the extrapolated 
impedance at zero frequency (Rq) which becomes real. From 
the Cole-Cole model the impedance at infinite frequency Roo 
is also extrapolated. (Roo) is also real and represents the total 
body water resistance. (Ri) defined as the ICW resistance can 
be calculated from (RoRoo )assuming (Ri) and (Re) are in 
parallel. Fig. 6. 
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Fig. 6. Equivalent circuit model for the extracellular, intracelluar and cell 
membrane spaces. 



B. Estimating fluid volumes from bioimpedance data 
The impedance of a cylindrical shaped object is: 






ui^ 



(3) 



Where: 

Z: impedance in (Ohm) 

L: length of the cylinder in (cm) 

V: Volume of the cylinder in (cm^) 

p : Resistivity in (ohm. cm) 

In bioimpedance measurements, L is assumed to be the 

separation distance between the detection electrodes and V is 

the volume to be estimated. However biological tissues are 

not homogenous; in fact they are considered as a suspension 

of conduction and non conducting materials. The apparent 

resistivity of suspensions can be calculated by Hanai mixture 



PROC. CAIRO INTERNATIONAL BIOMEDICAL ENGINEERING CONFERENCE 2006 < 



theory (Eq.4) Hanai (1968) describes the apparent resistivity 
(p) of a mixture containing a nonconductive material with 
volumetric concentration (C) in a conducting media with 
resistivity (Po) [8] 

P-—^ (4) 

(1-C)'/' 

At low frequency ECW is considered as the conduction 
material, and ICW which is surrounded by cell membrane is 
the non conducting material .The volumetric concentration in 
Eq.3 can be written as : 

V 

^ prw 



C,,=l- 



(5) 



Where Vtty is the total tissue volume. 

At high frequency, the injected current will penetrate both 

ECW and ICW spaces, therefore we can write: 



C, 
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V -\-V 
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V 
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(6) 



The impedance at zero frequency, Zq is extrapolated form 
Cole mathematical model. Zq is real and describes the 
resistance of the ECW i.e (Zo=Re) 

Using the basic impedance Eq.l and using the volumetric 
concentration at LF Eq.5, the volume of the ECW can be 
calculated by: 



F. 
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(7) 



Using the extrapolated values of Ro, Roo and the volumetric 
concentration at HF Eq.6, the volume of the ICW 
compartment can be calculated by: 

' ^ ^0 r. .. F. 
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Where 



K,= 



(9) 



and Pecw, Picw are the resistivity of the ECW,ICW in 
(Ohm. cm) respectively which can be found in literature. 

II. Methodology 
A. Bioimpedance sensor 

The main objective of this work is to develop a non-invasive 
method based on local multifrequency bioimpedance 
measurements that allow us to determine the fluid distribution 
and variations during haemodialysis. For this purpose a 
measurement system was developed by our research group 
[9]. Fig.7 




Fig. 7. Bioimpedance sensor developed at our research lab 
(UPC/Spain). 



The main specifications of this bioimpedance sensor are listed in 
Tabid. 

Table .i 

SPECIFICATIONS OF THIS BIOIMPEDANCE SENSOR 



Characteristic 


Value 


Measurement range 


of module 


5(^)<|Z|< 200(a) ±2.0% 


Measurement range 


of phase 


0.5(deg) < 1^1 < 90(deg) ± 2.0% 


Electrical current 




650 laA 


Frequency range 




lOkHz-lMHz 


Power type 




Battery operated (low power) 
Li-Ion 


Operating time 




8 hr at SR=(1 meas / 4 min, 6 
freq.) 


Weight 




Lightweight : 140(g), including 
cables, connectors and electrodes 


Dimensions 




60 X 60 X 25(mm) 


Cables 




20(cm) coaxial 


Electrodes 




Surface electrodes (SM Red Dot 
(TM) Ag/AgCl 2560) 


Measurement type 




Local Tissue Impedance 
spectroscopy 



B. Experimental protocol 

Clinical measurements were done in the Nephrology 
department/Medical center Pare Tauli /Sabadell/Spain. In 
these clinical measurements 60 haemodialysis sessions were 
included with average duration of 4 hours for each HD. 
9patients (3M,6F) participated in these clinical trails (6-8 HD 
session for each patient). The basic patient data are shown in 
TableJI. 

Table, ii 

PATIENT DATA 



Patient 


Age (year) 


Sex 


Height (cm) 


Dry Weight (kg) 


1 


33 


F 


165 


60 


2 


46 


F 


159 


67 


3 


77 


F 


151 


66.5 


4 


45 


F 


157 


50 


5 


55 


M 


182 


73 


6 


51 


F 


162 


59.5 


7 


81 


M 


165 


81 


8 


46 


F 


152 


56 


9 


53 


M 


173 


90 


10 


42 


F 


151 


64.5 



For each patient, three impedance sensors were attached to 
the arm, abdomen, and leg. At each segment the separation 
between the voltage detection electrodes was 10cm, and the 
injection electrodes were placed 5cm away from the voltage 
electrodes. Local tissue impedance module and phase were 
continuously recorded with 6 frequencies 

(10,20,50,100,220,500 kHz) with SR of 1 meas/4min. 
Continuous measurements of the total blood volume change 
(ABV) and heamatocrit (Hct) were obtained by an optical 
HCT analyzer (Crit-Line System). Intra dialysis clinical 
parameters such as: arterial blood pressure, ultrafiltration 
rate, and ultrafiltration accumulated were recorded every 20 
minutes. 

III. Results 
A. Pre and Post-HD states 

Fig. 8 shows the impedance module recorded from the leg at 6 
frequencies (10,20,50,100,220,500kHz) during the HD 
treatment session. To compare the variations in the 
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impedance between the beginning and the end of the HD 
session, the impedance module and phase are plotted in the 
frequency domain for the pre and post-HD states as shown in 
Fig.9 and 10. At both states the measured impedances are 
fitted into the Cole-Cole model and plotted in the complex 
plane in order to extrapolate the impedance values at zero and 
infinite frequencies Fig. 1 1 




1 00 1 50 

Time (min) 



Fig. 8. Bioimpedance Module recorded in the HD time. (Leg). 




Fig.9. Bioimpedance Module .Pre-HD and post-HD in the frequency domain. 
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Fig. 10. Bioimpedance Phase .Pre-HD and post-HD in the frequency domain 
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Total tissue volume (TTV) is estimated by measuring the 
circumferences and the separation between the detection 
electrodes. TTV decreases during HD the reduction 
percentage of the TTV is assumed to be equal to the total 
reduction of the body volume produced by the total 
ultrafiltration volume. For the pre and post-HD states the 
ECW and ICW volumes are estimated using the extrapolated 
values of Ro and Roo as explained in Eq.7 and 8. Total tissue 
water (TTW) is defined as the summation of ECW and ICW 
volumes. The average values of the 6 HD sessions of each 
patient at pre and post-HD are listed in Table.III and IV. 

TABLE. Ill 

ESTIMATED VOLUMES FOR THE PRE AND POST-HD IN (cm^) 

Patient TTV(O) TTV(f) ECW(O) ECW(f) ICW(O) ICW(f) TTW(O) TTW(f) 



641.32 150.79 118.24 304.28 349.55 

837.76 177.89 147.77 432.32 447.60 

574.96 125.35 103.70 205.75 193.34 

449.17 114.19 93.99 187.08 186.11 

722.31 121.04 93.56 190.72 179.14 

586.78 555.24 115.53 86.84 163.44 161.59 

643.24 125.11 111.14 246.01 242.59 

466.86 434.54 102.44 80.65 233.68 225.51 

477.38 459.81 97.87 80.21 157.51 165.62 



TABLE. IV 

ESTIMATED VOLUMES FOR THE PRE AND POST-HD IN (cm^) 

ECW(0) ECW(f) ICW(0) ICW(f) ECW(0) ECW(f) ICW(0) ICW(f) 

"tTWCO) "xTWCf) "tTWCO) °TTW(f) " TTV " TTV " TTV " TTV 



1 


685.63 


2 


857.28 


3 


599.81 


4 


468.73 


5 


752.73 


6 


586.78 


7 


656.71 



454.97 


426.26 


610.22 


576.30 


330.99 


273.14 


301.28 


261.39 


311.72 


243.48 


278.98 


218.61 


371.12 


340.55 


336.12 


276.09 


255.38 


228.93 



Patient 



1 34.08 30.22 65.94 80.21 22.20 

2 29.43 25.94 70.57 77.40 20.78 

3 37.82 37.95 62.21 70.82 21.52 

4 37.88 35.96 62.12 71.19 24.40 

5 38.81 38.43 61.21 73.54 16.08 

6 41.49 39.81 58.51 73.94 20.17 

7 34.20 32.95 65.80 70.98 19.14 

8 30.47 29.29 69.53 81.65 22.01 

9 38.28 35.25 61.72 72.10 20.43 



18.73 


43.54 


52.80 


6.20 


17.69 


50.22 


53.22 


2.17 


18.46 


35.25 


34.70 


3.98 


20.96 


40.14 


41.55 


3.94 


12.95 


25.33 


24.79 


3.85 


16.01 


28.57 


29.89 


5.19 


17.33 


37.26 


37.53 


1.97 


18.55 


50.20 


52.16 


6.58 


17.38 


33.08 


36.27 


3.50 



B. Continuous estimation of ECW and ICW during HD 

For continuous monitoring of the variations in ECW and 

ICW during HD session, the same procedure explained 

previously is repeated at each point of the measurements 

instead of the beginning and end state of HD only. 

Fig. 12 shows the impedance in the complex plane and the 

Cole-Cole fitting for all the points from pre to post HD. 
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Fig. 12. Continuous representation of the Bioimpedance in the complex. 
Measured points (dots) and the fitting Cole-Cole model (solid line). 



Fig. 1 1 . Bioimpedance in the complex plane at the pre-HD and post-HD. 
Measured points (dots) and the fitting Cole-Cole model (solid line). 
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Relative changes in the ECW and ICW are calculated by: 

Q,,(o/o) = 100 ^^^(^)-^^^(^ = ^) (10) 

'''''^ ECW(t = Q) ^ ^ 



C,c^(%) = 100 



ICW(t)-ICW(t = 0) 



ICW{t = 0) 
Fluid balance is defined as 

ECW 
And the relative change in FB is calculated by : 
ICW , , ICW , ^, 



'ICW 
ECW 



ECW 



ECW 



ICW 



(11) 



(12) 



(13) 



{t = 0) 



ECW 

LEG— Fluid Change(%) 




Fig. 13. Relative changes (%) for the ECW, ICW and fluid balance during the 
HD session. 

IV. DISCUSSION 
Bioimpedance spectroscopy (BIS) is used to evaluate the 
volumes of the ECW and ICW compartments. This approach 
includes the injection of electrical current normally below 
l(mA) at a range of frequencies. The questions that emerge 
here are: 'what is this frequency range? How low and high 
frequency is needed? and how many frequency points ?' 
The answer to these questions is : In order to evaluate the 
volumes of ECW and ICW we need to inject the electrical 
current in a range that goes as low as few (Hz) and as high as 
few M(Hz) with as many very high frequency samples as 
possible. However it is clear that this is impossible due to 
technical and health safety reasons especially in clinical 
applications. Instead, bioimpedance is measured at a selected 
frequency range and number of points and fitted into a 
mathematical model in order to extrapolate the impedance 
values at zero and infinite frequencies. 

In this work the bioimpedance spectroscopy (BIS) method is 
proposed to estimate ECW and ICW volumes using a self 
designed BIS sensor during HD treatment sessions using 6 
frequencies (10,20,50,100,220,500 kHz). Results show a 
significant decrease in the measured impedance during HD 
induced by water removal via ultrafiltration Fig. 8. The Cole- 
Cole model was able to fit the measurement points at the 6 
selected frequencies with high accuracy as shown in Fig.9-1 1 
and the values of Rq and Roo were extrapolated with high 
confidence. 



The total relative changes in the estimated total tissue water 
(TTW) and intracellular water (ECW) showed a significant 
decrease between the pre and post-HD for all patients as 
shown in Table. Ill and Fig. 14. 
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Fig. 14. Total relative changes (%) for the ECWa n d TTW after the HD 
session. 

The estimated intracellular volume (ICW) did not show a 
constant decreasing behavior as recorded in ECW and TTW 
as shown in Table.III and Fig. 15. 



% change in ICW volume 



(%) 




Fig. 15. Total relative changes (%) for the ICW after the HD session. 

The ECW ratio to the total tissue volume (ECW/TTV) and 
the ECW ratio to the total tissue water (ECW/TTW) were 
evaluated for the pre and post-HD states. The same was done 
for the ICW. A clear decrease in the ECW ratios were always 
recorded between the pre and post-HD states, however an 
increase was mostly recorded for the ICW ratios as shown in 
Table.IV and Fig. 16-17. This can be explained by the fact 
that ultrafiltration involves water removal which mostly 
comes from ECW space pace. 
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Fig. 16. Total relative changes in ECW/TTW and ECW/TTV after the HD 
session. 
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Relative change in ICW/TTW and ICW/TTV 
■ ICW/TTW 



□ ICW/TTV 
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Fig. 17. Total relative changes in ICW/TTW and ICW/TTV after the HD 
session. 

V. CONCLUSIONS 

Bioimpedance spectroscopy method was used to estimate 
the ECW, ICW and TTW volumes for patients during 
haemodialysis treatment sessions using 6 frequencies. The 
Cole-Cole mathematical model was used to fit the measured 
points in order to extrapolate the impedance values at zero 
and infinite frequencies. Precise matching was recorded 
between the measured points and the fitting model which 
indicate an accurate impedance sensor in the measurement 
frequency range. ECW and TTW volumes and the ratios 
ECW/TTV ECW/TTW showed a clear decrease in response 
to ultrafiltration during HD in all patients. However the 
behavior of ICW volume variation was not the same in all 
patients. The ratios ICW/TTV and ICW/TTW showed an 
increase between the pre and post-HD states which indicates 
that water is removed more from the extracellular space 
during ultrafiltration. 
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